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SCATTERING  FROM  A  VEGETATION  LAYER  WITH  AN 
IRREGULAR  VEGETATION  SOIL  BOUNDARY 


INTRODUCTION 

This  research  report  presents  a  theory  for  analyzing  the  nature 
PURPOSE  of  radar  wave  scattering  from  certain  types  of  vegetation.  The 
vegetation  is  simulated  by  a  continuous  random  medium,  and  use 
is  made  of  a  first-order  renormalization  technique  to  calculate  the  radar  backscatter 
coefficient.  The  influence  of  an  irregular  vegetation  soil  interface  has  also  been  con¬ 
sidered,  using  a  noncoherent  approach. 


In  a  previous  report  ,1  a  derivation  was  presented  of  the  radar  back- 
BACKGROUND  scatter  coefficient  from  a  half  space  of  random  medium  using  a 
first-order  renormalization  solution  for  the  scattered  wave  and  an 
isotropic  correlation  function  for  the  random  dielectric  fluctuations.  Recently,  Fung 
solved  the  problem  of  scattering  from  a  vegetation  layer  by  using  a  scalar  first-order 
renormalization  approach.2  In  his  solution,  however,  he  did  not  consider  the  existance 
of  a  rough  vegetation  soil  boundary.  He  did  consider  an  anisotropic  correlation  func¬ 
tion  in  which  the  horizontal  variation  is  different  from  the  vertical. 


Tsang  and  Kong  solved  the  problem  of  volume  scattering  from  a  half  space  random 
medium  that  contains  lateral  and  vertical  fluctuations.3  A  radiativ  transfer  approach 
was  used  to  calculate  the  backscattering  cross  sections  up  to  secor. .  order  in  approx¬ 
imation.  This  enabled  the  cross  polarized  terms  to  be  obtained. 


There  are  two  important  practical  applications  for  developing  and  analyzing  various 
radar  scattering  theories.  The  first  application  is  radar  image  simulation  of  terrain 
features.  In  this  problem,  the  radar  system  parameters  and  terrain  parameters  are 
known  and  used  to  calculate  a  radar  response  in  the  form  of  a  gray  tone  or  density. 


*R.A.  Hevenor,  Backscattering  cf  Radar  \.jves  By  Vegetated  Terrain,  U.S.  Aimy  Engineer  Topographic  Labor*- 
iorie*.  Fort  Belvoir,  VA.  ETL-0105,  June  1977,  AD-A047  669. 

2 

A.K.  Fung,  “Scattering  From  a  Vegetation  Layer,"  IEEE  Trantactiont  on  Geoscience  Electronics,  Vol.  GE-17, 
No.  1,  January  1979. 

3L.  Ttang  ana  J.A.  Kong,  “Radiative  Transfer  Theory  for  Active  Remote  Sensing  of  Half  Space  Random  Media,” 
Radio  Science,  VoL  13,  No.  5,  September-0 ctober  1978. 
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The  scattering  theories  can  be  used  to  compute  the  radar  backscatter  coefficient, 
which  in  turn  is  used  to  calculate  gray  tone.  Using  scattering  theories  in  this  type  of 
application  is  straightforward,  even  though  a  solution  for  any  one  particular  scattering 
problem  may  be  extremely  complicated. 


The  second  application  is  in  the  field  of  remote  sensing  of  terrain  in  which  the  sensor 
responses  must  be  used  to  determine  various  terrain  parameters  Using  scattering 
theories  for  this  application  is  not  straightforward. 


However,  there  are  two  important  uses  of  scattering  theories  that  bear  directly  on 
remote  sensing.  The  first  use  is  a  parameter  sensitivity  study.  The  theory  can  be  used  to 
analyze  the  influence  of  various  vegetation,  terrain,  and  radar  parameters  upon  the 
sensor  response.  Such  parameters  as  surface  roughness,  soil  moisture,  vegetation  height, 
and  density  could  be  varied  one  at  a  time  to  determine  the  influence  on  the  sensor 
response.  This  type  of  analysis  should  lead  to  determining  what  radar  parameters  are 
most  sensitive  to  certain  terrain  parameter  changes.  This  type  of  analysis  assumes  the 
existance  of  scattering  theories  that  have  been  developed  and  compared  with  existing 
experimental  data. 


The  second  use  is  to  analyze  the  radar  response  for  two  different  types  of  terrain 
features  to  see  if  the  two  features  could  be  distinguished  from  each  other  on  an  image. 
Once  again,  this  would  assume  the  existance  of  scattering  theories  that  have  been 
developed  and  tested  against  experimental  data.  These  applications  provide  the  in¬ 
centive  for  developing,  analyzing,  and  testing  various  scattering  theories. 


In  this  report,  the  geometry  of  the  scattering  problem  to  be  solved  and  the  basic 
technique  used  for  the  solution  will  be  discussed.  In  the  analysis  section,  the  derivation 
of  the  necessary  equations  will  be  provided.  In  the  results  section,  the  resulting  theory 
will  be  compared  with  existing  experimental  data,  and  a  study  on  the  sensitivity  of  the 
input  parameters  will  be  provided. 


In  this  report,  the  rationalized  MKS  system  of  units  is  used.  A  line  under  a  symbol 
will  be  used  to  represent  a  vector  quantity.  A  double  line  under  a  symbol  will  be  used 
to  represent  a  dyadic.  A  list  of  the  most  important  symbols  is  provided  at  the  end  of 
this  report. 
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In  figure  1 ,  the  scattering  geometry  of  the  vegetation  problem  it  shown.  A  plane  wave 
with  a  time  harmonic  of  exp  (j<ut)  is  incident  from  free  space  at  an  angle  onto  a 
layer  of  vegetation.  The  mean  thickness  of  the  vegetation  is  L.  The  vegetation  soil 
boundary  is  considered  to  be  randomly  rough  according  to  the  tangent  plane  approx¬ 
imation.  The  vegetation  is  simulated  by  a  continuous  random  medium  in  which  e  (.1 ) 
and  o  (JL  )  represent  the  three-dimensional  random  dielectric  and  conductivity  fluc¬ 
tuations,  respectively.  These  fluctuations  consist  of  the  sum  of  an  average  and  a  fluc¬ 
tuating  component.  The  standard  deviations  of  the  fluctuations  arr  represented  by  rjj 
and  t?2.  The  angle  of  refraction  of  the  mean  wave  in  the  random  medium  is  9(. 


Z 


FIGURE  1.  Scattering  Geometry. 
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The  soil  below  the  vegetation  represented  by  medium  3  is  assumed  homogeneous 
with  a  complex  propagation  constant  ka.  The  magnetic  permeability  for  all  three 
media  is  assumed  to  be  that  of  free  space.  The  electric  field  (  )  incident  onto  the 

vegetation  layer  can  be  written  as  follows: 

E,  -  jajig  +  a2ay  +  a3SJe*jko<xsinMcos0,) 

where  §x,  *y,  and  »,  are  unit  vectors  in  the  x,  y,  and  z  directions,  respectively. 
The  constants  at,  a2,  and  a3  are  arbitrary,  allowing  for  the  consideration  of  both 
horizontal  and  vertical  polarizations.  A  first-order  renormalization  method  will  be  used 
to  calculate  the  mean  and  scattered  waves  in  the  random  medium.  A  solution  will  be 
developed  first  for  the  case  where  the  vegetation  soil  boundary  is  a  plane  interface. 
The  irregular  boundary  will  be  considered  in  a  noncoherent  manner  afterwards.  The 
dielectric  and  conductivity  fluctuation  terms  (  e'  ( I )  and  o' (!)  are  considered  as 
being  generated  by  statistically  homogeneous  random  processes.  The  means  and 
correlation  functions  of  the  random  processes  are  defined  as  follows: 

<*'(!)>  =  <o’(r)>  =  o 


<e'(I)e'(l')>  =  <o'(J)o'(  l')> 


c-  lx  —  x*  1/8*  g* ly  “  y'l/8y  e-l*  -  t'l/fi* 


where  £x,  £  ,  and  £s  are  the  correlation  distances  in  the  x,  y,  and  z  directions, 
respectively.  The  correlation  functions  have  been  chosen  to  be  anisotropic.  This  re¬ 
presentation,  with  unequal  correlation  distances,  is  believed  to  be  closer  to  reality  than 
an  isotropic  correlation  function.  This  is  because  the  size  of  vegetation  scatterers  in 
a  horizontal  plane  is  not  the  same  as  the  size  of  the  scatterers  in  a  vertical  plane.  The 
mean  wave  in  the  random  medium  is  determined  from  the  bilocal  approximation  of 
the  Dyson’s  equation: 


IVxVx  -  kjl  <£(!)>  -  <IHU')>  <E(l')>  •  J?  ( I,  If )  sir’ 0 


where  <  £  (X)  >  is  the  mean  wave  in  the  random  medium. 


{(!)-  -jWMo^o'd)  +  c^Mo^e'd). 

k!  ■  +  w2iiQt0*k. 

r  d,  I* )  is  the  dytdic  Green's  function. 

V'  is  the  volume  of  the  random  medium. 


In  the  next  section,  plane  wave  solutions  will  be  sought  to  the  Dyson  equation  using  an 
infinite  space  dyadic  Green's  function.  It  should  be  noted  that  the  mean  wave  is 
located  in  the  integrand,  making  any  solution  very  difficult.  Once  the  mean  wave  has 
been  calculated  and  the  appropriate  boundary  conditions  have  been  matched  for  the 
mean  waves  in  all  three  media,  then  the  scattered  wave  in  the  vegetation  layer  can  be 
calculated  from  the  following  equation: 


[Vx  Vx  -  kj]  E.  (I)  -  {  (I)  <1  (  r  )> 


where  E^  (  r  )  is  the  scattered  wave. 


The  mean  wave  acts  as  a  source  term  for  the  scattered  wave,  which  will  be  computed 
using  a  Fourier  transform  technique.  This  in  turn  will  enable  the  scattered  waves 
in  air  to  be  determined.  The  necessary  boundary  conditions  will  be  matched,  and  the 
backscatter  coefficient  will  be  calculated  for  horizontal  and  vertical  polarizations. 
The  influence  of  the  rough  boundary  between  the  vegetation  and  the  soil  will  be  con¬ 
sidered  apart  from  the  volume  scattering  solution,  using  the  tangent  plane  method.  The 
backscatter  coefficient  for  rough  surface  scattering  will  be  modified  by  the  attenuation 
through  the  vegetation.  This  result  will  then  be  added  to  the  volume  scattering  solution 
to  obtain  a  final  answer  for  the  backscatter  coefficient.  An  elementary  permittivity 
model  will  be  developed  that  will  relate  certain  parameters  of  the  random  media  to 
the  parameters  of  actual  vegetation.  A  discussion  of  results  section  will  follow  in  which 
the  theoretical  results  are  compared  with  actual  experimental  data.  Also,  a  parameter 
sensitivity  study  will  be  conducted  on  the  theory  to  determine  the  influence  of  various 
parameter  changes  upon  the  final  result. 
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ANALYSIS 


In  this  section,  the  necessary  mathematical  derivations  will  be  provided  to  enable 
a  scattering  model  for  vegetation  to  be  obtained.  First,  a  solution  for  the  mean  wave 
will  be  presented  and  then  the  scattered  wave  will  be  calculated. 


MEAN  WAVE  The  first  step  in  obtaining  a  solution  to  the  Dyson  equation 

SOLUTION  is  to  write  the  dyadic  Green’s  function 


G0  (!,!') 


(1) 


where  i  is  the  unit  dyadic,  and  G0  (I,  l!  )  is  the  scalar  Green’s  function  that  satisfies 
the  following  equation: 

(V2  +  k2  )G0  (L,l! )  =  6(I-r')  (2) 

The  solution  for  equation  (2)  is  usually  expressed  in  terms  of  R  =  ll  -  l'  L  How¬ 
ever,  this  particular  form  is  not  useful  when  working  with  an  anisotropic  correlation 
function  that  is  expressed  in  rectangular  coordinates.  We  shall  therefore  seek  a  solution 
to  (2)  that  uses  rectangular  coordinates.  Let  G0  (  I,  l!  )  take  the  following  form: 

G0(I,l')  =  (-~3  /  dk^fJt)  efc  (3) 


where 


k 


+  k  a 
y-y 


Substituting  (3)  into  (?)  results  in  a  solution  for  G0  (k  ). 


S0U) 


k2  —  k2  —  k2  _  k2 
■  Rx  *y 


(4) 
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When  (4)  is  placed  in  (3)  and  integration  is  performed  in  the  complex  kz  plane, 
G0  (!,JL  )  becomes 


G»ai,)  "  5r/*»  f.  ^ 

.  exp[itkx(x  —  x*) +  kv(y  —  y*)  —  k,  iz-z'  1)1 

K 

K  =  V  *1  -  kj  -  k ; 

Transforming  (5)  into  polar  form  results  in 
i  r  2ff 

G“(X'I')  =  s?  J  dk  /  d9 

O  •'o 

^  exp[j{k(x  -  x')cos0  +k(y-y')sin0  ->/k2-  k2  lz-z'l}] 

*  (6) 

k  =  yk;+kj  kx  =  kcos0  ky  =  ksinfl 

When  (6)  is  used  in  (1),  the  dyadic  Green’s  function  becomes 

2ir 

1(1,1')  =  -^-  f  dk  f  d0~^=^{C(k,0)B(k,0)/k2-i} 

o  o  V  a 

(7) 

•  exp[j  { k(x  -  x')cos0  +  k(y  -y')sin0 ->/ k2  -  k2  Iz-z'  l}] 

where 

C(k,  0  )  =  ax  (jkcos0)  +  ay  (jksin0)  +  azf2(k) 

_B(k,0)  =  ix(-jkcos0)  +  Sy(-jksin0)  +  3zfi(k) 
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and 


f,(k) 


f2(k) 


j  \Zk\-  k1 

when 

z  >  z' 

-l/f? 

when 

z  <  z' 

when 

z>  z' 

j  A2-*2 

when 

z  <  z' 

Plane  wave  solutions  to  the  Dyson  equation  take  the  following  form: 

<E(I)>  =  Ae-JSe  *-t 


The  vector  A  can  be  obtained  by  matching  boundary  conditions.  The  mean  wave 
is  seen  to  propagate  with  an  effective  propagation  constant  Jse,  which  must  be  deter¬ 
mined.  When  the  above  equation  is  placed  in  the  Dyson  equation  along  with  the  dyadic 
Green’s  function  given  by  (7),  the  following  equation  is  obtain  when  the  cross  corre¬ 
lation  terms  between  the  dielectric  and  conductivity  are  ignored: 


[  Vx  Vx  -  k^J]  Ae-'-'Sc  *l  =  o 


(8) 


£  =  —  I - -  M 

k2o  ~  (2tt)3  ~ 


M  ~  jn(v2V22  y*  dk  J 


(9) 


d0 


d(k)e~lx-x'  1  / fixe~ly  -y'  l/«ye-l*-*'  l/«. 


.  e+jk e  l  C(k,0)B(k,»)/kJ-I]  exp  [j  {k(x-x')cos0 

+  k(y  -  y ')sin0  -  y/  k*  -  k2  I  z  -  z'  I  |  J 


(10) 


d(k)  =  k  /i/k2  - 


The  integral  in  li  is  allowed  to  be  over  all  space.  For  the  case  where  £r  is  very  small 
this  should  be  a  good  approximation  except  for  the  points  extremely  close  to  either 
boundary.  By  carrying  out  the  integration  in  X',  one  obtains  a  result  for  My  that 
represents  the  ijill  element  of  the  dyadic  .M. 


Mij  =  4KP*CyC*  f  dk  /  <Wd(k){  g(k) 

•  [Ft  (k,  0)^,0)  +  ^]  +  h(k)  [F2l(k,e)F2J(k,  0) 

+  kjfiy]}/  |l+£j(keR+kcose)2l 

•  [l+£y(key+ksin0)2]|  (11) 


where 


g(k)  =  - 

+7^)  h(k) 

DZ  V  A 

-ykj-k2)-! 

Fjjfk,  0)  =  jkcos0 

F12(k,  0)  =  jksinfl 

F13(k,0)  =  j/k2-k2 

F21(k, 0)  =  jkcos0 

F22(k,0)  =  jksin0 

F23(k,0)  =  -j/k2-k2 

6.,  =  1  when  i  -  j 

=  0  when  i  =£j 

K  =  jir(T?2r?2  -k2r?J)/k2 
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The  form  of  the  incident  wave  dictates  that  ksy  =  0  and  kex  =  kosin0(.  This 
result  makes  the  following  elements  of  the  £  dielectric  tensor  become  equal  to  zero: 

*12  =  C2I  =  C23  C32  “  0 

The  above  result  is  easily  shown  by  considering  a  transformation  of  (1!)  back  to 
rectangular  coordinates  and  recognizing  that  tie  integrand  is  even  in  ky.  Carrying 
out  the  indicated  differentiation  in  equation  (8)  and  writing  the  result  in  matrix 
form  produces 


k2  -k2£ 

ez  oil 


-k  k  sinfl.  -  k2£ 

ez  o  i  o  1 3 


k2  +  k2sin0.-k2e,, 

ez  o  i  o  22 


k  k  sin0.  -  k2?. ,  0 

ez  o  i  o  1 3 


kjsin2e,-k^„ 


-ike 


A  e‘Ji® 

X 


A  e*j*e  - 


A  e-jS*  ** 

Z 


In  forming  the  ebove  matrix,  use  has  been  made  of  the  fact  that  key  -  0  and 
kex  =  k0sin6j.  Now,  only  solutions  for  kez  are  needed.  Two  solutions  can  be 
developed  for  kei,  one  for  a  horizontally  polarized  wave  and  one  for  a  vertically 
polarized  wave.  For  a  horizontally  polarized  wave,  one  has  Ax  =  Ax  =  0  and 
Ay  =A  0.  The  following  equation  can  be  used  to  determine  kez  for  this  case: 


|k2  +  k2sir20.  -  k2?,,}  A  e'jk®  ’x  *  0 

l  ez  o  1  o  22’  y 


Since  the  term  outside  the  brackets  is  not  zero,  this  means  that  the  quantity  inside  the 
brackets  must  be  zero. 


k  =  ±k  v/ sin20,  - 
ez  o  v  i  22 


The  above  result  is  not  an  explicit  solution  for  kez  since  this  quantity  also  appears 
in  the  integral  of  M2i.  A  first  approximation  for  kez  can  be  obtained  by  letting 
*?i  *  t?2  =  0.  This  represents  the  case  where  there  are  no  random  fluctuations  at 
all. 


k^>=  VkJ-kJsin2®, 


The  above  value  can  be  used  to  compute  M22>  which  in  turn  can  be  used  to  cal¬ 
culate  a  new  value  of  k„  that  will  be  called  kh.  The  minus  sign  on  the  square  root 
is  chosen  over  the  plus  sign  to  consider  waves  propagating  in  the  minus  z  direction. 


When  computing  kh  in  (15),  the  expression  for  k^  is  used  to  calculate  M22. 
It  is  interesting  to  see  that  even  when  k,  is  real,  kh  st  'l  comes  out  complex  so  that 
the  mean  wave  decays  as  it  propagates  into  the  medium.  This  decay  nas  been  explained 
as  resulting  from  multiple  scattering.  It  is  not  clear,  however,  how  much  multiple 
scattering  is  being  considered.  For  a  vertical  polarized  wave,  Ax  0,  Ay  *  0, 
and  Az  ^  0.  This  leads  to  the  following  determinant: 


Once  again  a  first  approximation  for  kei  can  be  obtained  for  the  case  where  the 
random  dielectric  and  conductivity  fluctuations  disp pear  (tj,  =  tj2  =  0). 

k<°>  =  _ 


This  value  for  kez  can  be  used  to  calculate  the  elements  of  the  dielectric  tensor. 
These  elements  are  used  in  (16)  to  compute  a  new  value  for  k„,  which  will  be 
called  kv. 


The  sign  associated  with  the  square  root  has  been  chosen  as  minus  in  order  to  consider 
waves  propagating  in  the  minus  z  direction.  The  effective  propagation  constant  for 
the  mean  wave  has  been  determined,  and  now  the  amplitude  must  be  calculated 
by  matching  appropriate  boundary  conditions.  The  total  mean  electric  field  in  air 
can  be  written  as: 


E,(r)  -  j  Ca1Sx  +  a2Sy  +  a3az]  e,koBCO,®t 


+  [R,ax  +R2§y  +R3§i1  e**'ko*COifl‘  }  e-jk-*»,nfl* 

z  >  0  (17) 

The  first  bracketed  term  in  (17)  is  the  incident  wave,  and  the  second  bracketed  term 
is  the  reflected  wave.  The  unknowns  in  (17)  are  represented  by  Rj,  R2,  and  R3. 
However,  it  will  not  be  necessary  to  obtain  an  explicit  solution  for  them  since  they 
are  not  needed  in  determining  the  scattered  waves.  The  total  mean  electric  field  in  the 
vegetation  (E2(r))  can  be  written  in  the  following  form,  using  previous  results: 


l2(r)  =  {  T2ayeP1*e,q‘*  +<'T1»X  +  +  V2Sye'p2VJ,2* 


-L  <  z  <  0  (18) 

where 

Pi  +jqt  =  -jkh 

Pj  =  Im(kh)  q,  =  -Re(kh) 

also 

P2  +  jq2  =  -jky 

p2  =  Im(ky)  q2  =  -Re(ky) 
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In  (18),  both  upward  and  downward  waves  have  been  considered.  The  unknowns 
are  represented  by  the  amplitudes  Tj,  T2,  Tj,  Vj,  V2,  and  V3.  An  explicit 
solution  for  each  of  these  is  required  in  terms  of  propagation  constants,  medium 
characteristics,  and  layer  thickness.  The  explicit  solution  is  needed  to  compute  the 
scattered  waves.  An  expression  for  the  mean  wave  in  the  homogeneous  soil  medium 
can  be  written  as 


E3(I)  =  (W,ax +W2§y +W3aj^kox,ln0ieik3*co,®3 


z  <  -L  (19) 


In  the  soil,  the  mean  wave  propagates  in  the  minus  z  direction  and  the  constants 
Wj,  W2,  and  W3  are  the  unknown  amplitudes.  Once  again,  no  explicit  solution 
will  be  required  for  these  amplitudes  since  they  are  not  needed  to  compute  the 
scattered  waves.  To  compute  the  six  amplitudes  of  the  mean  wave  in  the  vegetation, 
the  following  boundary  conditions  are  used: 


E 

E 


lx 

"  E2, 

at  z 

ly 

=  E2y 

at  z 

0 

0 


JSl*  _ 

3Ely 

3E2i 

3E2y 

3y 

dz 

3y 

dz 

dE, 

1  X  _ 

3E1*  . 

dE2x  _ 

9E2z 

3z 

dx 

dz 

3x 

at 


=  0 


at  z  =  0 


E2x  =  E3x  at  2  =  ~L 
E2y  =  E3y  at  Z  =  ”E 


3E2i  3E2 

3y  dz 


9E3z 

3y 


-L 
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5,  .  Jh± 

dz  dx 


dE,  dE, 

_li  -  at  z  *  -  L 

dz  dx 


Dl*  =  D2i  at  z  =  0 


°2*  =  D3i  at  z  =  -L 


=  -o-i 

D2  =  evE2x-2x  +  ehE2ySy  +e¥E2iat 


where 


q2  -  p2 

— ^ — - — -  evaluated  at  0{ 


=  0° 


„2  2 

Q.  -  p. 

— - i -  evaluated  at  0. 

,2 .  1 


=  0° 


B,  = 


3-3 


where  e3  is  the  dielectric  constant  of  the  soil.  Two  divergence  conditions  will  also 
be  used,  along  with  the  above  boundary  conditions. 


V  •  I,  =0  and  V  •  I3  =  0 


When  the  equations  for  the  mean  fields  given  by  (17)  through  (19)  are  placed  in 
the  boundary  conditions,  the  result  is  12  equations  and  12  unknowns.  Explicit 
solutions  are  only  required  for  the  six  amplitudes  associated  with  the  mean  wave  in 
the  random  medium.  Solving  for  these  six  values  yields  the  following  results: 


2jk0a12a2cosflt 

a2la12-al  1  a2  2 


(20) 
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(21) 


~2jkna1a1  ^osflj 
a2lai2"®Ua22 


Ti 


•ikobi2  {  2a1cosfl1  +  sinfl^l  -c0/ev)[a3  -t-a^anfl)]} 
b21b12_bllb22 


(22) 


V 


1 


-jkpt),  i  {2a,cosfl1-t-sinei(l  -c0/fv)[a3  +a>  tangi)| 

b2 1 bl 2~bl 1 b2  2 


(23) 


t3 


V 


3 


_ 1_  |  e3kosin0, 

ev  Ie-P2Lp-Jq2L_eP2Leiq2L]  1“^^ 

•  (T. e'PlL(  iq2L  +V,epjLeiq2L)  -e  ePjLejq2L 

*  1  O 


•  [a3  -  tan0,  (T,  f  V,  -a,)l 


(24) 


e0/ev  [a3  -  tanfl,  (Tj  +  V,  -  a, ))  -  T,  (25) 


The  parameters  us’d  iu  the  abo^c  equations  are  defined  below: 
a,,  =  (jk3cos03  -pj  -jqj)e'PlLe'JqiL 

ai 2  3  Ok,cosfl3  +  n.  *jq1)ePlLejqiL 

a2 1  =  jkpeosOj  +  Pj  +jqt 

a2  2  =  jkpCCsflj  -  p,  —  jq3 

bn  '=  eP2Le'Jq2L(p2 +jq2 -a) 
b  1 2  =  -cP2Lejq2L(p2 +jq2 +a) 
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where 


b2 1  ■  P2  ♦Mj  +jk0tcoa#l  +  si*a,tandl(l-«0/«v)l 
b22  “  "  |Pj  +j*i2  -jk0lcoad,  ♦  stouten#, (l-«p /*,)l| 

*  =  j  {kjcostfj  +  k*sina0i(l-eJ/e¥)/(kJcoaflJ)} 


Now  that  the  mean  waves  have  been  fully  determined,  the  scattered  waves  can  be 
calculated.  The  scattered  waves  in  the  upper  medium  (air)  will  then  be  used  to  com¬ 
pute  the  backscatter  coefficient. 


SCATTERED  WAVE  In  the  random  medium,  the  scattered  or  incoherent 
SOLUTION  field  is  calculated  from  the  following  equation: 

[Vx  Vx  -  kJlE^U)  *  $ (I)  <§(!)>  (26) 

For  our  problem,  the  mean  wave  <  IL  (  X  )  >  is  given  by  £2(X)  as  shown  in  (18). 
A  superscript  2  is  used  with  J,(2)(!)  to  indicate  clearly  the  scattered  field  in  the 
random  medium.  A  solution  for  this  scattered  field  can  be  obtained  by  using  a  two- 
dimensional  Fourier  transform. 


-•1,(X)  *  w-  /  <27> 

where 

.k,  =  kx§,  +  ky8y  It  =  xflx  +  y*y 

£,<*,.*)  =  G./k,,*#,  +  G„<*,.«)l,  ♦  G(1(iSt,  z)i. 


When  the  complete  expression  for  the  mean  wave  in  the  random  medium  as  given  by 
(18)  is  placed  in  (26),  a  term  of  the  form  {  (  r )  exp(-jk0xsinfl,)  on  the  right  side 
results.  This  term  will  b3  written  as  a  two-dimensional  Fourier  transform. 
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1.1‘HK'III 


yiV".-. 

fey* 

r 


S(kt>  z)  =  y* S lft|U)e'ik°x,infl‘e‘J*» 

f(De'Jkox,,n0‘  =  — J—  /  dktS(^,  z)ej*‘ **»  (28) 


When  (27)  and  (28)  are  placed  in  (26)  and  the  result  is  put  in  matrix  form, 


"(kJ-kJ-D  l) 

-kxk(  jk,D, 

V 

~kxky 

(kJ-kJ-D?)  jk,D, 

G.y 

« 

fyU) 

Jk,D.  (k^kj-kj) 

f,(z) 

where  D,  represents  the  differential  operator  d/dz.  The  quantities  on  the  right 
side  of  (29)  are  defined  below: 

fx(z)  =  S(k1,z)(T1ep2*ei,2*  +  v,e*Pj*  e'Jq*1] 

fyU)  =  S(kt,z)lT2ep»*ejq»*  +  V2e*Pllejp»*] 

fz(  z  )  =  S(kt,  z)[T3ep2*  ejqjZ  +  V3e'p**  e'Jq2‘) 


Solutions  for  the  quantities  G,x ,  Giy ,  and  G„  can  be  obtained  by  solving  the  three 
differential  equations  in  (29)  using  the  method  of  variation  of  parameters. 
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'■.A-*-  [/  h(i^k;'di] 


G**(-t’  z>  = 


(k,A,  +  kvB,)  jk;t  (k.Ai+Mi) 
ki  ~  k! 


•  e1* 


»•  + 


s*;2*? 


J  £kxf(z}  -  kyh(z)j  ^‘‘ds 

&  /  M-H 

•'-L  L  J  (32) 


k'  *  >/k2  -  k2  -  k* 

f(z)  -  fx(*)lkj-kj) -kykxfy(z)+jkxDIfI(z) 

h(  z  )  -  ky  -k,-~ fxU)  -  (k,2  -  kj)fy(z)  - jkyDIfI(z) 
*x 

+  (ki  *kJ-kJ>f.(z) 


The  quantities  A),  AJt  Bj,  and  B}  are  not  functions  of  z  and  at  present,  are 
unknown.  Although  the  solutions  for  Gtx,  Gty,  and  G„  are  rather  formidable  in 
appearance,  it  will  be  found  after  some  mathematical  manipulations  that  a  solution 
will  emerge.  The  three  components  of  the  scattered  electric  field  in  the  upper  medium 
(z  >  0)  can  be  written  in  the  following  form: 
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(33) 


EV>>  ■  f 


E(1)(£)  ■ 

Iy  11  (2tr) 


where 


f  &xAyOLtyJ{k*x *  kyy  ‘k»*x) 


E(1)(n  c  JL,  /*dk  A  (jLW(k*x  *  k*y  “  kllI) 

« v  (2rY  J  1 


(34) 


(35) 


kl«  “ 


The  parameter  kt  t  w r  obtained  by  taking  any  one  of  the  field  components  given 
above  and  putting  it  i  the  free  space  scalar  wave  equation.  The  superscript  1  is 
used  to  refer  to  the  fiel<  the  upper  medium,  which  is  air.  Therefore,  >(-£-)  would 
indicate  the  x  component  of  the  scattered  electric  field  in  air.  The  components  of  the 


scattered  electric  field  in  the  soil  (z  < 

-  L)  can  be  written  as  follows: 

E<3>(_L)  = 

tx  v  7 

(2ir)2  J 

^*SiktCx(kt)e1(kxX  +  kry  +  k3*l) 

(36) 

E(3)(L)  = 

(S*)2  J 

dk  r.  (k  )e*(l  *x  +  kyy  *  k3*l) 
t  y 

(37) 

E(3)U)  = 

IZ  9 

_L  / 
(2*)2  J 

• 

dJc^C  (kt)ej<k*x  +  kyy  +  k3**) 

(38) 

where 

k3z  = 

v/k2-k 

2  -  k2 

X  Ky 

The  superscript  3  refers  to  the  computation  of  the  scattered  fields  in  the  soil.  The 
expression  for  k3l  is  obtained  by  putting  any  one  of  the  field  components  into  the 
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scalar  wave  equation,  which  has  a  propagation  constant  k3.  The  unknowns  associated 
with  the  scattered  waves  are  represented  by  Ax,  Ay,  Az,  A1(  A2,  Bt ,  B2,  Cx,  Cy, 
and  C2.  The  only  unknowns  for  which  an  explicit  solution  is  needed  are  Ax,  Ay> 
and  Az,  Since  the  only  interest  is  in  computing  the  backscattered  far  field  in  air, 
complete  solutions  for  the  scattered  fields  in  the  other  mediums  are  not  needed.  The 
boundary  conditions  that  the  scattered  waves  must  satisfy  are  provided  below: 


EU)  =  E(2)  at  z  =  0 

IX  IX 

E(1)  =  E(2)  at  z  =  0 
iy  »y 


3E(1)  3E(1)  3E(2)  3E(2) 

sz  _  »y  _  sz  _  «y 


3y 

3z 

3y 

3z 

3E(1) 

sx 

3E(1> 

SZ 

3E(2) 

SX 

3E(2) 

SZ 

3z 

3x 

3z 

3x 

e(2)  = 

IX 

e(3) 

11  IX  ' 

=  -L 

p(2)  - 

11  »y 

E(,3y>  at  z 

=  L 

3E(2) 

_ il—  _ 

3y 

3E(2) 

_ UL.  - 

3z 

3E(3) 

_ U— 

3y 

3E<3> 

_  _ a_ 

3z 

3E(2) 

IX 

3E(2) 

iz  - 

3E(3) 

IX 

3e£> 

3z 

3x 

Sz 

3x 

at  z  =  0 


at  z  =  0 


at  z  =  -  L 


at  z  =  -  L 


The  boundary  conditions  given  above,  along  vith  the  divergence  equations  in  the 
two  homogeneous  media  allow  ten  independent  equations  to  be  formulated. 


kxAx  +  kyAy  =  k,zAz 


(39) 
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-xAx^kosin0i>  0)  +  -y  Ay(kosin0j,  0) 


2  cos d{  (jk0)e'jk°R 
4ffR 


+  atAl(k0sinfli(  0) 


(49) 


Where  R  is  the  distance  from  the  origin  of  the  coordinate  system  to  the  field  point 
where  E,f  is  required.  It  can  now  be  seen  that  solutions  for  Ax,  Ay,  and  Az  are 
only  needed  for  values  of  kx  =  kosin0j  and  ky  -  0.  When  the  appropriate  expressions 
for  the  derivatives  are  substituted  into  equation*  (45)  through  (48)  and  we  let 
ky  =  0,  the  10  equations  (39)  through  (48)  become. 


k*Ax  = 

klzAz 

(50) 

= 

kaA 

(51) 

A,  - 

I, 

A,  +  A,  +  - \-r- 

1  2  2jk;k  l 

(52) 

n 

>. 

< 

B  +  B,  +  - r-r— 

1  2  2jk;k,2 

(S3) 

Cx  = 

[*■' :i  + 

(54) 

Cx  = 

eJk  3iL  +  Bjeik,L 

(55) 

JklzAy  = 

-«i'2  - 

(56) 

+  kxAz  = 

kJA.-kjAj  +  0+k;/k^)l4 

,  k?,A  A  ,  kf(0) 

T^(Ai"A2) - 

Kz  *z 

(57) 
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-k;B,eJk*L  +  k;B2ejk*L  =  k3lCye'jk3*L 
A2(k;3  +  k*k;)e-Jk‘L  -  A1(k1,,  +  ^WkiL  +kxfz(-L) 


(58) 


*  k:vjk3*L(k3zcx-kxcr) 


(59) 


The  quantities  Ip  I2,  I3,  and  I4  used  in  the  above  equations  are  integrals  in  z  and 
are  defined  below: 


The  functions  f(z),  h(z)  and  k  z  can  be  evaluated  at  ky  =  0,  before  integration 
takes  place.  When  all  the  above  equations  are  used  to  solve  for  Ax,  Ay,  and  Az, 
the  following  results  are  derived: 


Ax(kx,o) 


bjfz(-L)  +  b2f2(o)  + 


f(z)  |  b5eik ** 


+  b6e'^kz*  |  dz 


(60) 
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»,  *  2k1,(k,Jtk;kJ»J/I> 

b2  *  -ki*kxd+52)/a8 

b3  “  a7  ( 1  "*■ 

b4  =  a5/(2j?8k;kj) 

bS  =  b3  +  b4 
b6  b3  b4 

If  the  receiver  in  the  far  field  is  sensitive  to  a  unit  polarization  vector  £r(  then  the 
received  field  (ER)  will  be 


where  <-r  in  general  has  three  components  (§r  -  3xerx  +  3yeiy  +  3,,erz).  The  next 
step  in  calculating  the  backscatter  coefficient  is  to  determine  the  statistical  average 
of  ERER,  using  (49). 


<ErEI>  " 


4k2  cos20j 
16tt2R2 


cer*x  <AxAx*>  +  2Relerxery 


•  <AxAy*>  |  ♦  2R=(e,xe*  <A,A*>  I  +  ^ y 

•  <A,a;>  +  2ReleIye*(  <  Ay  A*  >  1  +  e„e* 

•  <  A, A*  >  |  (63) 


*  t 
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The  buckets  < - >  around  ER  ER  are  used  to  indicate  the  calculation  of  the 

statistical  average.  The  possibility  of  erx  and  efl  being  complex  is  anticipated, 
with  y  always  being  real.  What  is  required  now  is  the  computation  of  each  of  the 
six  terms  inside  the  brackets  of  (63).  The  determination  of  <AyA*>  will  be  con¬ 
sidered  first:  y 


<  AyA*  >  =  k2k*2  M0M: 


MoMo  J  dz  J  dz  ^  dx 

_•  "  •• 

J  dX  /  dy  j  dy'<H-t  )f*d'  )  > 


f(z,  z')  exp(-2jk0(x  -  x')sin0jj 


where 


<t(D  «*(!')>  =  (w2^+kjn}>r,x-x'  l/eVly**y'  “  z'  l'fiz 


f(z,z')  =  a,  a  *T2  X*  eD  1  zeD  t 1  +  ajjajvje0  izeD2*' 


+  a.T^e0!^'0^'  +  a,T2V 


*aDii  -D  j  z' 

,  e  1  e  1 


+  a,  VaSjTj  e°2*  eD>z'  +  a^VjVj  e°2*  e°2z' 
+  ®i  V2Tj  eD  2Z  e*D  2*'  +  a" ,  V2  V*  eD  2Z  e‘D  tz' 

+  af^TjT*  e'D 2*  eD  iz  +  aJX2V*e'D2Z  eD2*' 

+  T2 X*  e’D  2 1  e"D  2 z '  +  x  V*  e‘D  2 z  e'D  l 1 ' 

+  V2aJX*cDizeDiz'  +^V2V2*  eD«zeD2z' 


+  V2Xje-Di*eD2«'  +  V2V*  eD  j  *  eDTJ 
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where 


Dj  *  pt  +  j(q,  -  k|)  and  D2  ■  -<Pj  +  jqt  +  jk') 


°  ikJ  kltO+ai)  +  k^l-a^ 

The  form  of  f(z,  z')  given  above  appears  to  be  very  complicated.  However,  each  of  the 
16  terms  in  f(z,  z')  consist  of  simple  exponentials  in  z  and  z'  and  therefore  can  be 
integrated  easily.  It  should  be  remember  when  computing  D}  and  D2  that  pt 
and  qt  are  real,  but  k't  will  be  complex.  Making  the  substitution  that  u  ■  x  -  x' 
and  v  =  y  -  y'  and  transforming  the  x  and  y  integrals  into  integrals  in  u  and  v 
produces 


<  AyA*  > 


o  ••  mm  rn 


dy' 


f(z,  z’)  e-2jkou,lne‘  e'lu 


I / fix  - 1  v  I  /fiy  -  Is  -  *'  l/C* 
c  c 

(65) 


The  integrals  in  x'  and  y'  appear  somewhat  meaningless.  These  integrals  actually 
represent  the  illuminated  area  in  the  xy  plane,  since  it  is  physically  unrealistic  to  have 
backscattered  energy  from  a  portion  of  the  surface  that  is  not  illuminated.  Considering 
the  integrals  in  x'  and  y’  to  form  the  illuminated  area  (Aj)  and  carrying  out  the 
integrations  in  u  and  v  will  yield  the  following  result: 


<  AyA*  > 


m„m; 

(i  +  4k2cWe1) 


f  dz  f  dz' 

-L  -L 


f(z,  z')e 


- 1  *  -  *'  l/£* 


(66) 
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Consider  now  a  typical  term  of  f(z,  z')  which  is  of  the  form  Aealeb>  where  A, a, 
and  b  are  not  functions  of  z  or  z’,  and  make  a  transformation  of  variables  from 
z  and  z'  to  n,  ■  z  -  z'  and  z"  3  z’.  Then,  the  results  of  carrying  out  the  integration 
for  this  one  term  becomes 


f  dz  f  dz' 

•'-L  -L 


Ae  e 


6b*'e- It  - »'  l/£i  _ 


Afi. 


(a  +  b) 


1  -8tb  +  g^a+b*44**  1/6x)-(H-Cta)e'L(‘  +  b) 
(1  +fi,a)  (1  -C,b) 


l-8,a  +  8I(a  +  b)e‘L<1/e*  +  b)-e'L(,+  b)(l+b8I) 
(l  +  8rb)  (1  -  «xa) 


When  the  answer  for  the  integration  in  z  and  z'  given  above  is  used  for  each  of  the 
16  terms  in  ?(z,  z'),  then  the  final  result  for  <  AyA*  >  can  be  written. 


<  Ay  A*  > 


4gx«y  W<o2My2  4 

(1  +  4kjfi2  sin20.) 


16 


Y]  (c„ 


+  d„) 


1  *czdn  +ez<cn  +d„)e'L<c"  +  lliz)  -  (1  +8tcn)e~L(Cn  *  d»> 


(i+filCn) 


n  *  1 


(67) 


l-g*cn  +  gr(cn +dn)e'L(dn  +  1/fil)-e  L(c"  +  dn)  (1 +dnet) 
(1  +  g*dn)  (1  -  £rcn) 


The  values  for  the  A„'s,  the  cn's,  and  dn's  are  provided  in  appendix  A  and 
simply  come  from  the  expression  for  ?(z,  z').  Using  the  methodology  for  computing 
<  AyAj  >,  one  can  caluclate  all  the  remaining  terms  in  (63).  All  of  these  other 
terms  are  given  in  appendix  A.  An  expression  for  the  radar  backscatter  coefficient 
(a°)  can  be  written  in  terms  of  <  ER  ER  > 
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The  subscript  v  on  o°  is  used  to  indicate  a  volume  scattering  result  from  a  plane 
layer  of  random  medir..  No  consideration  for  rough  surface  scattering  is  given  in  a?. 
Using  (63)  and  the  expression  for  the  incident  wave  given  previously,  one  can  write 
a  final  result  for  a° : 

ao  =  kocos2ei 

ff(aia*  +  a2a*  +  a3a*) 

•  erx  ®xx  +  2Re  [erx  ery  a  yl 

+  2Re  [efx  e*  «XJ  +  er2y  «yy  +  2kosin0, 

•  Re  Ky  e«  °S«y  /  kl  *  1  +  e  x*  er*.  J  <69> 

where 


“xx  = 

<\a*>/a, 

«xy  = 

<AXA*>/A, 

«xz  = 

<aX>/a, 

«yy  = 

<  Ay  A*  >  /  A, 

* 

<  A*A*>  /  A, 

Consider  now  the  form  of  o°  for  horizontal  and  vertical  polarizations.  The  following 
parameters  are  used  to  describe  a  wave  that  is  transmitted  with  horizontal  polarization 
and  horizontal  polarization  is  received. 


0 


f  y 


e  =0 

vr* 


For  these  parameters,  the  backscatter  coefficient  can  be  given  the  additional  sub¬ 
scripts  of  HH  to  indicate  horizontal  polarization  transmit,  and  horizontal  polarization 
receive. 


HHv 


kj  ayy  cos2  0,/  it 


((70) 


The  case  of  vertical  polarization  transmit,  vertical  polarization  receive  can  be  char¬ 
acterized  as  follows: 


at  =  cos0j 


rx 


COS0j 


a 


2 


0 


ry 


*  0 


a3  =  stnfij 


et,  =  sin0, 


The  backscatter  coefficient  associated  with  these  parameters  can  be  given  the  addi¬ 
tional  subscripts  VV. 


k2  cos20j 


'  VVv 


jaxx  cos20j  +  2Re  [axl  cosOjSindj] 


+  sin2®i 


(71) 


If  a  result  is  computed  for  the  cross-polarized  backscatter  coefficient  (HV  or  VH), 
the  term  will  disappear.  The  reason  for  this  is  that  the  particular  elements  of  the  dyadic 
M,  which  would  yield  cross-polarized  terms  in  the  mean  wave,  are  all  zero.  Next, 
let’s  consider  the  influence  of  an  irregular  vegetation-soil  boundary. 
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In  this  section,  we  will  consider  what  must 
be  done  to  equations  (70)  and  (71)  to  include 
the  influence  of  a  rough  ground  surface.  It  is 
expected  that  the  influence  of  the  rough  surface 
would  be  greater  when  the  angle  of  incidence 
is  small.  Abo,  as  the  vegetation  height  or  density 
gets  larger,  less  scattering  is  expected  from  the 
ground  surface  below.  In  what  follows,  the  horizontal  and  vertical  polarizations  will  be 
considered  separately. 


MODIFYING  THE  VOLUME 
SCATTERING  RESULTS  TO  IN¬ 
CORPORATE  THE  INFLUENCE 
OF  AN  IRREGULAR  VEGETA¬ 
TION  -  SOIL  BOUNDARY 


Consider  a  horizontally  polarized  wave  incident  from  free  space  onto  a  layer  of  vegeta¬ 
tion  that  has  an  average  thickness  L.  The  interface  between  the  vegetation  and  soil 
will  be  considered  as  randomly  rough  in  such  a  way  that  the  tangent  plane  approxi¬ 
mation  is  applicable.  The  radar  backscatter  coefficient  (o£  H )  will  be  considered  as 
the  sum  of  a  term  resulting  from  surface  scattering  and  a  term  resulting  from  volume 
scattering. 


<£«  "  0^HSexPI"4tt.iLsec^.i1  +  °HHt  (72) 

In  equation  (72),  o^HS  represents  the  backscatter  coefficient  fora  randomly  rough 
surface  with  a  guassian  distribution  of  surface  heights.  The  subscript  s  indicates  sur¬ 
face  scattering.  The  quantity  H  «  is  multiplied  by  a  decaying  exponential  in  which 
ote  {  is  the  imaginary  part  of  the  effective  propagation  constant  and  x  is  the  true 
angle  of  refraction  for  the  mean  wave  in  the  vegetation.  The  second  subscript  1  on 
a,  |  and  \pt  i  is  used  to  indicate  horizontal  polarization  since  these  parameters  will 
have  different  values  for  vertical  polarization.  The  effective  propagation  constant 
(k# , )  can  be  obtained  from  kh . 

k*  *  k3  +  k3  =  k  3  +  k3sin30. 
tl  cz  ex  hoi 

If  we  now  let  ke,  =  -  jatsI  and  in  place  of  kh  we  put  jp,  -q,,  then  we  can 

solve  for  0e  j  and  ae  j . 


0,1  *  P,i  cos(*el/2) 
a,l  *  P*ii  sin(tf>e,/2) 


i 


i 

I 


i 


»• 
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Where  pt{  and  0(1  are  defined  below: 


J  4pJqJ  +  (qf  +  kJsin29,-p})2J  * 


♦el 


tan 


-t 


2ptqj 


q?  +  kjiina«  -p} 


The  planes  of  constant  phase  for  the  mean  wave  in  the  random  medium  are  used  to 
calculate  an  expression  for  sec^e  t . 


see*,, 


✓5  sfn20i  +  pj(gelcos71  -ot^sin^) 
P,  <0etcosY,  -ft.jSinT,  ) 


Where  p,  and  are  given  below: 


Pi 


4ajbjsin4  6{  +  [1  -(a*  -bj)sin2  0J2 


Vl* 


-  Vi  tan 


-i 


2a,btsin2 


1  -  (aj  -  bj)sin2fl, 


•l- 


o  "el 


& +  «Ji 


b»  - 


ko“el 


fill +  «e2l 


Many  derived  expressions  are  available  for  the  backscatter  coefficient  from  a  randomly 
rough  surface  using  the  tangent  plane  method.  The  following  equation  will  be  used:5 


^IHS 


g2  | 

T~r  -  1  -Sin2  *el  12(1  'goCOs2  *el> 

4m*  cos  | 

exp  [-tan2  ^el/(4m2)l  (73) 


-  sin2  t  (1  +  g2)J 


1 


5R.  A,  Hevenor,  Backscattering  of  Electromagnetic  Waves  From  a  Surface  Composed  of  Two  Types 
of  Surface  Roughness,  U.S.  Army  Engineer  Topographic  Laboratories  Fort  Belvoir,  Virginia. 
FTL -TR-71-4,  October  1971,  AD-737  675. 
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.ihikiC  ■' 


The  quantity  R0  in  (73)  represents  the  Fresnel  reflection  coefficient  and  can  be 
computed  as  follows: 

0.1  cos*.i  s*"2  ~K\ 

0.1  cos  *.i  +  sin2  *tl 


In  calculating  R0 ,  we  have  neglected  the  effect  of  the  imaginary  ports  of  the  pro¬ 
pagation  constants.  The  term  m(  represents  the  ratio  of  the  standard  deviation  of  the 
surface  fluctuations  to  the  correlation  distance.  The  quantity  g0  is  defined  by  the 
following  expression: 


8o“  1- 


20.1  cm 


v'0?  -  0J,  “"J 


A  final  equation  can  now  be  written  for  oH  H  that  considers  both  the  volume-scat¬ 
tering  and  surface-scattering  effects. 


4m2  cos4  i^ol 


1  -  sin2  i//eJ  12(1  -g0cos2 


-  sin2  ^el  (1  +  g2)I  J  exp  (-tan2  ^el/(4mJ)] 
.  exp[-4ael  Lsec  ^el  j  +  k2  ayy  cos 2di/ir 


(74) 


In  the  same  manner,  a  complete  solution  for  vertical  polarization  can  be  obtained: 


°v  v 


4m2  cos4  \lit2 


ro  sin2  cos 

1  +  - - -  (2r0sin^eJ 

To 


•  K  , 

+  ro  cos  ^e2  ]  +  - - sin  \j/t2  cos3  \pt2 
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•  exp  [-tan2  \pt2/ (4m2))  exp  [-4  «e2  L  sec  \J/t2 ) 


k2  cos20j 


oxx  cos20j  +  2Re  [cosflj  sin0j  axz)  +  azz  sin20. 


The  new  quantities  introduced  into  the  above  equation  are  defined  below: 


0e2  =  COS  (^t2^ 


at2  =  P* 2  sin  (0e2/2) 


The  quantities  pe2  and  0e2  are  given  as  follows: 


P*  2  = 


4P2  d2  +  (d2  +  ko  sin20i  ”  pj>2  |  W 


«e2  =  tan* 


^2 _ 

q2  +  k2  sin20t  -  p2 


k2  sin20j  +  p2  (P2t2  cos72  -  ae2  sin^)" 
p2  (/3e2  cos  72  -  ae2  smr2) 


sec  tf/e2  = 


4a2 b2  sin40.  +  [1  -  (a'  -  b')  sin'0.) 


2  l2\  „;_2n  i  2 


7 2  =  Zi  tan’ 


-i  r  2a2b2  sin2gi 


1  -  (a2  —  bj  )  sina01 


where  a2  and  b2  are 


Mei 


Hi  cos  ^e2  -  j?e2  sin2 

fi\  cos  *e2  +  0e2  v/^3  -  4  sin4  *e2 

{#»  Kl  +  /3e2  y/^^Tt2 } 


*2  " 


koae2 

/32,  +  a2, 

Ke2  e2 


Equations  (74)  and  (75)  are  the  final  results  for  the  radar  backscatter  coeffieient 
for  horizontal  and  vertical  polarizations.  Before  the  results  of  computing  equations 
(74)  and  (75)  are  shown,  an  elementary  vegetation  permittivity  model  must  be 
developed  that  relates  some  of  the  model  input  parameters  to  the  con  olex  dielectric 
constants  of  vegetation  and  water. 


To  determine  the  influence  of  various 
DEVELOPMENT  OF  A  VEGETATION  vegetation  parameters  (such  as  moisture 
PERMITTIVITY  MODEL  content)  upon  the  calculation  of  the 

backscatter  coefficient,  one  must  relate 
some  of  the  permittivity  parameters  in 
the  scattering  model  to  the  physical  parameters  of  the  vegetation.  Peake  and  Oliver’s 
model6  will  be  used  to  calculate  the  relative  complex  dielectric  constant  of  vegetation 
( *v  ): 


ev  =  (  F/2  )Re[  ew  ]  +  j(  F/3  )Im[  ew  ] 


^W.H.  Peake  and  TX.  Oliver,  The  Response  of  Terrestrial  Surfaces  at  Microwave  Frequencies,  Technical  Report 
AFAL-TR-70-301,  The  Ohio  State  University,  Electroscience  Laboratory,  AD-884  i06. 
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where  F  is  the  fraction  of  water  by  weight  in  the  vegetation;  Re[  2W  ]  and  Im[  ew  ] 
are  the  real  and  imaginary  parts  of  the  relative  complex  dielectric  constant  of  water 
(  ew  ),  which  can  be  written  as 


e  =  5  +  - — 

1  +  j(  1.8S/X  ) 

where  X  is  the  wavelength  in  centimeters.  For  particular  values  of  X  and  F,  we  can 
now  compute  ev.  With  a  knowledge  of  ey,  one  can  estimate  the  average  relative 
complex  dielectric  constant  (  )  by  using  the  following 


(V 


v  ev 


+  eA  VA 


)/v. 


Re[  ] 


where  Vy  is  the  volume  occupied  by  the  vegetation;  VA  is  the  volume  occupied 
by  air;  VT  is  the  total  volume  equal  to  Vy  +  VA ;  eA  is  the  relative  dielectric 
constant  of  air,  assumed  equal  to  1 .  The  variances  and  r?|  can  be  computed 
by  using  the  following  formulas: 


V 


2 

l 


V 


2 

2 


V,  (<  -  E.)2  +  VA  ( 


V,(c,  -  O,)2 


+  v*-; 


where 


av  =  -- weQF  Im[  ew  ]/3  and  e'  =  Re(ev] 

The  symbol  Rv  shall  be  used  to  designate  the  volume  ratio  Vy  /  VT . 

The  developed  model  for  the  radar  backscatter  coefficient  is  now  complete  and  cal¬ 
culations  can  be  made.  In  the  next  section,  computed  results  will  be  shown,  and  the 
theory  will  be  compared  with  some  existing  experimental  data. 
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DISCUSSION  OF  RESULTS 


In  this  section,  some  numerical  calculations  will  be  shown  for  the  theory  derived  in 
the  previous  section,  and  a  study  will  be  presented  of  the  influence  of  the  verious 
input  parameters  on  the  backscatter  coefficient.  Two  computer  programs  were  de¬ 
veloped  for  solving  equations  (74)  and  (75).  One  program  solves  for  equation  (74) 
and  the  second  solves  for  (75).  The  solutions  to  the  half-space  and  plane  layer  prob¬ 
lems  are  also  generated  for  comparison.  A  listing  of  the  computer  program  for  solving 
equation  (74)  is  given  in  appendix  B.  The  10  input  parameters  to  the  programs  are 

1 .  Fraction  of  water  by  weight  in  the  vegetation  (  F  ) 

2.  Volume  of  vegetation  divided  by  the  total  volume  (  Rv  ) 

3.  Correlation  distance  in  the  x  direction  (  2X  ) 

4.  Correlation  distance  in  the  y  direction  (  2y  ) 

5.  Correlation  distance  in  the  z  direction  (  2,,  ) 

6.  Mean  thickness  of  the  vegetation  layer  (  L  ) 

7.  Relative  dielectric  constant  of  the  soil  below  the 
vegetation  ( eg  ) 

8.  Conductivity  of  the  soil  below  the  vegetation  (  o3  ) 

9.  Frequency  ( f  ) 

10.  Ratio  of  the  standard  deviation  of  the  rough  surface 
fluctuations  to  the  correlation  distance  of  the 
fluctuations  (  m, ). 

The  output  of  the  computer  programs  is  the  backscatter  coefficient  in  decibels  as  a 
function  of  incidence  angle.  The  backscatter  coefficient  in  decibels  is  related  to  the 
backscatter  coefficient  as  follows: 

o°  (in  decibels)  =  101og10a° 

The  backscatter  coefficient  on  the  right  side  of  the  above  equation  is  computed  by 
(74)  or  (75).  The  following  discussion  centers  on  figures  2  through  30,  which  show 
the  results  of  computing  equations  (74)  and  (75).  Figures  2  and  3  come  from 
Ulaby  and  Bush  and  provide  pertinent  ground  truth  data  associated  with  the  experi¬ 
mental  measurements.  Figure  4  comes  from  Cihlar  and  Ulaby7  and  provides  a  re¬ 
lationship  between  soil  moisture  and  relative  complex  dielectric  constant.  Figures  5 
through  15  provide  a  comparison  of  the  developed  theory  with  experimental  data 
taken  from  a  cornfield  by  Ulaby  and  Bush.8 


7 

F.T.  Ulaby  and  T.F.  Bush,  Com  Growth  as  Monltores  by  Radar,  The  University  of  Kansas  Center  for  Research, 
Inc.  RSL  Technical  Report  117-57,  November  1975, 

o 

F.T.  Ulaby  and  J.  Cihlar,  Dielectric  Properties  of  Soils  as  a  Function  of  Moisture  Content,  The  University  of 
Kansas  Center  for  Research,  Inc.,  RSL  Technical  Report  177-47,  November  1974. 


FIGURE  2.  Com  Ground  Truth,  1974. 
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SOIL  MOISTURE  (GRAMS  PER  CM3) 


F.T.  Utaby  and  J.  Cihlar,  Dielectric  Properties  of  Soiit  u  «  Function  of  Moisture  Content,  The 
Univenity  of  Kinau  Center  for  Research,  lac.,  RSL  Technical  Report  177-47,  November  1974. 


Figure  4.  Representative  Dielectric  Constant  Values  as  a 
Function  of  Volumetric  Water  Content. 


The  soil  moisture  is  obtained  from  figure  2  for  a  particular  set  of  measurements  per¬ 
formed  on  a  given  date.  This  soil  moisture  is  used  along  with  the  curves  of  figure  4  to 
determine  the  relative  dielectric  constant  and  the  conductivity  of  the  soil.  Throughout 
all  comparisons  of  theory  with  experiment,  it  has  been  assumed  that  the  soil  type  is 
a  loam.  In  comparing  theory  with  experiment,  remember  that  certain  input  parameters 
to  the  theoretical  model  were  not  known  and  had  to  be  estimated;  whereas,  other 
input  parameters  were  known  from  the  ground  truth  data  collected  during  the  ex¬ 
periment.  The  unknown  input  parameters  are  Rv,  fix,  fiy,  fiz  and  m5. 

In  figure  5,  the  theory  is  matched  to  the  experimental  data  for  corn  that  is  at  a  height 
of  30  centimeters.  The  large  rise  that  occurs  in  o°  as  d{  goes  from  10°  to  0°  is 
indicative  of  a  rough  surface  effect.  In  this  case,  the  rough  surface  is  quasi-specular 
since  ms  is  given  such  a  small  value.  It  can  be  seen  that  to  match  the  theory  with  the 
experimental  data,  it  was  necessary  to  let  £x  be  different  from  fiy  and  to  let  fiz  be 
much  smaller  than  fix  or  fiy.  The  fact  that  fix  is  different  from  fiy  shows  an 
anisotropic  effect  in  the  horizontal  plane,  which  probably  arises  from  the  com  being 
planted  in  rows. 


Figure  6  matches  the  theory  to  experimental  data  for  com  that  is  2.3  meters  high.  It 
can  be  seen  that  a  good  match  is  obtained  for  8X  equal  to  By,  indicating  that  the 
anistropic  effect  in  the  horizontal  plane  has  essentially  disappeared  for  8.6  GHz.  The 
values  of  the  parameters  used  for  Rv,  8X,  £y,  8Z,  and  ms  in  figure  6  are  also  used 
in  figures  7  through  10  to  determine  whether  the  model  could  provide  a  correct 
prediction  of  a°  for  different  values  of  F,  soil  moisture,  and  vegetation  height. 

Figures  '  hrough  9  show  an  excellent  agreement  between  theory  and  experiment. 
Figure  lo  shows  an  excellent  agreement  between  theory  and  experiment  for  angles  of 
incidence  equal  to  and  greater  than  30°.  For  angles  of  incidence  less  than  30°,  the 
agreement  is  poor.  A  possible  reason  for  this  poor  agreement  may  be  due  to  the  rainfall 
that  came  prior  to  the  August  15  th  measurements.  The  rainfall  could  have  disturbed 
the  soil  surface  in  both  a  physical  and  an  electrical  manner  such  that  its  scattering 
behavior  is  no  longer  predictable  from  prior  values. 

Figures  1 1  and  12  show  an  attempt  to  match  the  theory  to  the  experimental  data 
for  frequencies  of  GHz  and  13  GHz.  It  can  be  seen  that  to  obtain  a  good  match, 

the  v  ...  v..;  fiy  fiy  must  be  altered  from  the  values  used  at  8.6  GHz.  This 

seems  to  indicate  that  as  the  frequency  goes  higher,  the  vegetation  medium  becomes 
more  complicated  and  the  anisotropic  behavior  becomes  more  pronounced. 
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Figure  6.  Comparison  of  Theory  with  Experimental  Data. 


Horizontal 


Figure  7.  Comparison  of  Theory  with  Experimental  Data. 


Polarization:  Horizontal 


Figure  8.  Comparison  of  Theory  with  Experimental  Data. 


Polarization:  Horizontal 
F  =  0.529  f  =  8.6  GHz 
R.  =  0.0002 


Figure  10.  Comparison  of  Theory  with  Experimental  Data. 


Polarization:  Horizontal 
F  =  0.839  f  = 


Comparison  of  Theory  with  Experimental  Data. 


In  figures  13  through  15  a  match  is  shown  of  the  theory  with  experimental  data  for 
vertical  polarization.  For  all  three  curves,  8X  must  be  made  unequal  to  2y  to  obtain  a 
good  match.  In  figure  16,  the  variations  are  shown  in  the  experimental  measurements 
of  o°,  which  can  occur  throughout  the  spring  and  summer  for  alfalfa.  Large  variations 
in  the  backscatter  coefficient  occur  prior  to  and  after  harvesting. 

Figure  1 7  presents  a  study  of  sensor-look  direction  with  respect  to  vegetation  planted 
in  rows.  The  two  parameters  and  oy  represent  the  backscatter  coefficient  when 
the  look  direction  is  perpendicular  and  parallel  to  the  rows,  respectively.  We  see  that 
when  0 ,  equals  zero  degrees,  o±  and  oy  are  equal.  However,  for  0.  greater  than 
zero  degrees,  is  greater  than  ojj.  The  theoretical  results  agree  only  partically  with 
the  experimental  results  given  in  figure  18,  which  comes  from  Batlivala  and  Ulaby.9 


Figure  19  provides  a  study  of  backscatter  coefficient  versus  layer  thickness  for  two 
angles  of  incidence.  For  the  0;  equal  to  zero  curve,  the  solution  for  o°  with  a  rough 
layer  differs  from  the  half-space  solution  by  approximately  16dB  (decibels)  for  a 
layer  thickness  of  0.5  meters.  As  the  layer  thickness  is  increased,  the  solution  for  o° 
at  0j  equal  to  zero,  approaches  the  half-space  solution.  The  0 j  equal  to  20°  curve 
also  approaches  the  half-space  solution  when  the  layer  thickness  is  increased.  In  this 
case,  o°  differs  from  the  half-space  solution  by  only  about  3.5dB  when  the  layer 
thickness  is  0.5  meters. 

Figure  20  presents  a  study  of  the  skin  depth  of  the  mean  wave  versus  incidence  angle 
for  three  different  frequencies.  For  a  horizontally  polarized  wave,  the  skin  depth  is 
taken  to  be  the  reciprocal  of  pj ,  which  was  derived  earlier  as  part  of  the  solution  to 
the  Dyson’s  equation.  It  should  be  remembered  that  the  mean  wave  decays  for  two 
reasons,  absorption  and  scattering.  For  a  frequency  of  8.6  GHz,  the  skin  depth  goes 
from  approximately  5  meters  at  0j  equal  to  0°  down  to  1  meter  at  0t  equal  to 
80°.  When  the  frequency  is  increased,  the  overall  level  of  the  curve  is  lowered  con¬ 
siderably,  but  it  does  not  drop  off  as  fast  with  increasing  incidence  angle. 

In  figure  21,  the  solutions  are  compared  for  the  half  space,  the  plane  layer,  and  the 
layer  with  a  rough  surface.  For  angles  of  incidence  between  0°  and  30° ,  the  rough 
interface  at  the  vegetation  soil  boundary  can  have  a  dramatic  effect  on  o° .  In  this 
case,  it  is  clearly  not  sufficient  to  use  a  plane  layer  model. 


9 

PJP.  Batlivala  and  P.T.  Ulaby,  The  Effect  of  Look  Direction  on  the  Radar  Return  From  a  Row  Crop,  National 
Aeronautics  and  Space  Administration,  Lyndon  B.  Johnson  Space  Center,  Houston,  Texas,  RSL  Technical  Report 
264-3,  May  1975. 
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14.  Comparison  of  Theory  with  Experimental  Data. 


Polarization:  Horizontal 


16.  Study  of  the  Experimental  Variations  of  the  o°  for  Alfalfa. 


Polarization:  Horizontal 


Figure  17.  Study  of  the  Sensor  Look  Direction. 
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SOURCE  :  PJ>.  Batlivala  and  F.T.  Ulaby,  The  Effect  of  Look 
Direction  on  the  Radm  Return  From  a  Row  Crop, 
National  Aeronautic*  and  Space  Administration, 
Lyndon  B.  Johnson  Space  Center,  Houston,  Texas 
RSL  Technical  Report  264-3,  May  1975. 


Figure  18.  Scattering  Coefficient  a°  as  a 
Function  of  Incidence  Angle  at  (a)  2.75GHz, 
(b)  5.25GHz,  and  (c)  7.25GKz.  Data  set 

*  1,  July  16.  1974. 
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19.  Study  of  the  Variation  of  a0  with  Layer  Thickness. 


Polarization:  Horizontal 
F  =  0.839 
R„  =  0.0002 


Study  of  the  Skin  Depth  of  the  Mean  Wave 
Versus  Incidence  Angle. 


Polarization:  Horizontal 


In  figures  22  through  30,  a  sensitivity-of-parameters  study  is  presented.  In  this 
study,  the  input  parameters  are  varied  individually  to  determine  the  overall  effect  on 
the  scattering  coefficient.  In  figure  22,  a  study  is  provided  of  o°  variations  with  F 
(the  vegetation  moisture  content).  When  0i  is  greater  than  20°,  the  larger  values  of 
F  produce  higher  levels  of  o°.  The  shape  of  the  o°  versus  6{  curve  does  not  change 
much,  but  the  overall  level  is  significantly  different.  At  approximately  0(  *  10°, 
a  crossover  of  the  curves  exists  such  that  the  curve  with  the  highest  moisture  content 
now  yields  the  lowest  value  of  o°.  The  reason  for  the  crossover  is  that  two  entirely 
different  mechanisms  are  responsible  for  scattering.  For  greater  than  20°,  the 
volume-scattering  mechanism  dominates  so  that  higher  moisture  in  the  vegetation 
results  in  larger  backscatter  values.  However,  for  angles  of  incidence  6t  less  than  10° 
when  the  mechanism  for  scattering  is  dominated  by  the  rough  surface  under  the 
vegetation,  the  higher  moisture  values  result  in  lower  o°  values.  This  is  because  the 
higher  moisture  values  provide  more  attenuation  of  the  mean  wave,  which  means  that 
less  is  available  for  scattering  from  the  rough  surface. 

Figure  23  presents  a  study  of  o°  variations  with  the  parameter  Ry .  The  smaller 
value  of  Rv  yields  a  larger  o°  value  for  small  angles  of  incidence.  This  is  again 
because  rough  surface  scattering  dominates  for  small  angles,  and  smaller  values  of 
Rv  mean  that  more  energy  gets  down  to  the  surface.  A  crossover  occurs  at  approx¬ 
imately  0j  =  15°  where  the  volume  scattering  result  begins  to  dominate.  Another 
crossover  occurs  between  dl  -  40°  and  0i  =  50°  such  that  for  angles  larger  than  50°, 
o°  falls  off  faster  for  the  larger  value  of  Ry .  The  reason  for  this  second  crossover  is 
possibly  because  that  for  the  larger  value  of  Rv ,  the  lower  interface  between  the 
vegetation  and  soil  no  longer  provides  a  contribution  for  backscattering. 

Figure  24  presents  a  study  of  a°  variations  with  L  (the  mean  thickness  of  the 
vegetation  layer).  For  small  angles  of  incidence,  the  smaller  value  of  L  yields  larger 
values  of  a°.  The  rough  surface  below  the  vegetation  is  dominating  the  return  and 
the  smaller  value  of  L  provides  a  lower  attenuation,  thus  making  more  energy  avail¬ 
able  for  scattering  from  the  surface.  A  crossover  point  occurs  around  0-  =  12°,  which 
indicates  that  volume  scattering  is  now  beginning  to  dominate  and  so  the  thicker 
layer  will  yield  a  larger  value  of  <r°.  Another  crossover  point  occurs  at  approximately 
0j  =  67°.  This  crossover  point  possibly  indicates  that  for  L  =  2  meters,  the  lower 
interface  is  having  no  influence,  but  for  L  =  0.5  meters,  the  lower  interface  still 
provides  a  contribution. 
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Polarization:  Horizontal 
Rv  =  0.0002  f  =  8  GHz 
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22.  Study  of  o°  Variations  with  F. 


Polarization:  Horizontal 


Figure  24.  Study  of  cr°  Variations  with  L. 


Figure  25  shows  a  study  of  o°  variations  with  soil  moisture.  The  soil  moisture  in¬ 
fluences  both  the  dielectric  constant  and  the  conductivity.  The  effect  is  shown  of 
doubling  the  soil  moisture  underneath  a  1 -meter  layer  of  vegetation.  Thus,  doubling 
the  soil  moisture  content  increases  the  level  of  o°  slightly  over  all  angles  of  incidence. 


In  figure  26,  a  study  is  presented  of  o°  variations  with  frequency  for  two  angles  of 
incidence.  We  see  only  a  slight  frequency  dependence  at  =  10°,  because  the  rough 
surface  scattering  is  independent  of  frequency  except  for  the  attenuation  portion. 
At  0j  =  30°  where  volume  scattering  is  playing  a  more  important  role, we  see  a  very 
significant  frequency  dependence. 

In  figure  27,  o°  variations  with  the  parameter  m,  are  shown.  This  parameter  re¬ 
presents  the  ratio  of  the  standard  deviation  of  the  rough  surface  undulations  to  the 
correlation  distance.  The  most  specular  surface  (m,  =  0.04)  yields  the  highest  value 
of  o°  at  0j  =  0°.  However,  for  this  surface,  a°  falls  off  very  fast  with  increasing 
incidence  angle  so  that  at  ei  =  20°,  the  rough  surface  effect  has  disappeared.  When 
m,  is  allowed  to  increase,  the  value  of  o°  at  =  0°  decreases.  Also,  as  m,  increases, 
the  rough  surface  influences  a°  over  a  larger  range  of  incidence  angles. 

In  figure  28,  o°  variations  with  £x  are  shown.  Notice  first  that  a  change  in  £x 
yields  virtually  no  influence  upon  0°  for  angles  of  incidence  less  than  10°.  For 
angles  of  incidence  greater  than  10°,  the  curve  associated  with  the  larger  value  of 
£x  is  higher  until  a  crossover  point  is  reached  around  6i  =  48°.  For  angles  of  in¬ 
cidence  greater  than  48°,  the  curve  associated  with  the  larger  value  of  fix  falls  off 
much  faster  than  the  curve  associated  with  the  smaller  value  of  £x.  Increasing  the 
value  of  £x  will  then  move  the  curve  upward  for  angles  of  incidence  less  than  about 
50°,  but  will  lower  the  curve  for  angles  of  incidence  greater  than  about  50°. 


In  figure  29,  o°  variations  with  £y  are  shown.  Once  again,  notice  that  a  change 
in  £y  has  virtually  no  influence  on  0°  for  angles  less  than  10°.  For  angles  of  in¬ 
cidence  greater  than  or  equal  to  20°,  an  increase  in  £y  results  in  an  increase  in 
0°,  Therefore,  increasing  £y  simply  increases  the  level  of  the  curve  for  angles  equal 
to  and  greater  than  20°. 


In  figure  30,  a°  variations  with  £z  are  shown.  For  angles  of  incidence  less  than  10°, 
changes  in  £z  have  no  influence  on  a°  owing  to  the  dominance  of  the  rough  surface. 
For  angles  of  incidence  greater  than  20°,  increasing  the  value  of  £z  simply  increases 
the  overall  level  of  the  curve  without  changing  the  shape.  It  can  be  seen  that  the  0° 
curve  is  sensitive  to  slight  changes  in  £t.  A  change  in  £z  of  only  a  fraction  of  a 
millimeter  produces  a  significant  change  in  o°.  This  sensitivity  may  make  any  attempt 
to  determine  £z  in  a  rigorous  experimental  manner  very  difficult. 
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Figure  25.  Study  of  o°  Variations  with  Soil  Moisture. 


Polarization:  Horizontal 
F  =  0.08  m,  =  0.04 

Rv  =  0.0002 


Study  of  o°  Variations  with  Frequency. 


Figure  27.  Study  of  o°  Variations  with  m,. 


Polarization:  Horizontal 
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29.  Study  of  a0  Variations  with  £y . 
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Figure  30.  Study  of  o°  Variations  with  £z . 


This  section  concludes  with  a  brief  discussion  on  the  limitations  and  difficulties  en¬ 
countered  in  developing  the  theory  presented  in  this  report.  It  appears  valid  to  simulate 
a  region  of  vegetation  with  a  continuous  random  medium,  although  it  is  not  certain 
as  to  how  well  the  first-order  renormalization  technique  does  in  solving  the  problem. 
It  is  not  clear  how  much  multiple  scattering  is  being  considered,  and  it  is  not  even 
clear  as  to  how  much  multiple  scattering  must  be  considered.  A  free  space  dyadic 
Green's  function  was  used  in  solving  the  Dyson’s  equation;  however,  what  should 
have  been  used  was  a  Green’s  function  applicable  to  a  layered  problem.  Such  Green's 
functions  are  very  complicated  to  develop  and  work  with.  Also,  it  is  not  clear  how 
much  the  final  result  will  change  if  a  more  complicated  Green’s  function  is  used. 


An  equation  for  the  radar  backscatter  coefficient  o°  was  developed  by  first  obtaining 
plane  wave  solutions  to  the  Dyson’s  equation.  Expressions  were  found  for  the  z 
component  of  the  effective  propagation  constant  for  both  horizontal  and  vertical 
polarizations.  The  mean  wave  was  then  used  to  calculate  the  scattered  wave,  which 
in  turn  was  used  to  compute  o°.  The  final  result  for  o°  indicated  the  necessity  to 
develop  a  permittivity  model  that  would  relate  some  of  the  permittivity  parameters 
in  the  scattering  model  to  the  physical  parameters  of  the  vegetation.  This  was  accom¬ 
plished  with  only  limited  success  since  a  rather  elementary  permittivity  model  was 
used.  The  final  result  for  o°  still  contained  five  input  parameters  that  were  unknown. 
The  five  parameters  are  2X,  £y,  fif,  Rv,  and  ms. 


Further  work  in  this  area  should  attempt  to  determine  correlation  functions  and 
distances,  as  well  as  Rv  and  m,  to  relate  clearly  theory  and  experiment.  An  aniso¬ 
tropic  correlation  function  was  used;  however,  this  did  not  result  in  a  depolarization 
term.  The  existance  of  depolarization  is  clearly  evident  from  the  experimental  data. 
The  reason  for  this  depolarization  is,  as  yet,  unknown.  A  depolarization  term  could 
be  obtained  by  computing  the  scattered  field  to  a  second-order  approximation.  It 
could  also  be  obtained  by  initially  allowing  for  an  anisotropic  random  medium.  At 
this  time,  it  is  unclear  which  approach  is  correct. 


'  •• 


CONCLUSIONS 

For  certain  types  of  vegetation,  such  as  corn,  the  irregular  vegetation  soil  bound* 
ary  dominates  the  backscattering  results  for  angles  of  incidence  between  0°  and 
20°.  Any  remote  sensing  of  surface  phenomena  beneath  vegetation  should  be 
done  in  this  angular  range. 


The  effect  of  the  rough  surface  boundary  between  the  vegetation  and  soil 
increases  with  decreases  in  frequency,  vegetation  moisture  content,  vegetation 
volume,  and  layer  thickness. 


Increasing  the  soil  moisture  content  increases  the  level  of  the  o°  curve  slightly. 


Using  different  correlation  distances  in  x  and  y  does  not  completely  explain 
the  effect  of  look  direction  on  scattering  from  a  row  cop. 


The  a°  curve  versus  incidence  angle  curve  is  sensitive  to  very  slight  changes 
in  the  correlation  distance  in  z. 


The  predictability  of  the  o°  curve  is  dependent  on  meteorological  phenomena, 
such  as  rain.  ' 


The  correlation  distances  and  the  vegetation  volume  ratio  do  not  stay  constant 
throughout  the  entire  growth  cycle  of  the  com  crop.  However,  once  the  crop 
matured,  these  parameters  remained  fairly  constant. 


Because  different  correlation  distances  were  required  to  match  the  experimental 
data  for  horizontal  and  vertical  polarizations,  a  more  correct  model  for  the 
vegetation  may  be  an  anisotropic  random  media  model. 


The  predictability  of  the  o°  versus  incidence  angle  curve  depends  upon  a  very 
detailed  knowledge  of  the  dielectric  fluctuations  of  the  vegetation  and  the  sur¬ 
face  roughness  properties  of  the  soil  below.  However,  such  knowledge  for 
particular  vegetation  features  does  not  exist  at  the  present  time.  This  detailed 
understanding  should  be  obtained  if  theoretical  models  are  to  have  ultimate 
usefulness  in  predicting  scattering  from  vegetation  features. 
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APPENDIX  A.  Definition  of  Terms  Involved  in  Computing 


<AyAj>,  <AXA;>,  <  A,A^>,  <AxAs*>,  <AxA^> 


Repeating  equation  (67)  produces 


.  _  4Cx8yeiA,(w2Mo^  +  kX)k.2kI2M0Mo 

y  y  =  (1  +  4kj£j  sin2®,) 

^  (  l-Mn+MCn^d„)e-L(cn^1/Ei)  -  (i+g,cn)e-L(cn  +  dn) 

2^  (cn+d„)  (  (»+fi.cn)  (1-M«) 


+ 


1  -  Ctcn  +  Cz(cn  +  dn)e"L(<‘n  *  ,/E*>--e'L(c"  +  dn)(l  +d„K 
_ _  (I  +w  (l-8iCn) 


The  values  for  An  ,  cn,  and  dn  are  defined  below: 


A> 

= 

cl 

s 

D. 

a2 

= 

a,T2a;v; 

C2 

= 

D« 

A3 

= 

w; 

C3 

= 

D, 

A4 

= 

C4 

= 

D! 

A5 

= 

«.v;n 

C5 

= 

°2 

\ 

= 

^1»1V2V2 

C6 

- 

D2 

A7 

= 

^2t*2 

C7 

= 

D2 
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Before  giving  the  values  for  A  n,  0n,  and  7n,  the  following  parameten  are  defined: 

D1  *  P2+jq2-jk; 

Da  *  “(p2 +Jq2 +jk') 

k,  -  b,  IVk.’-^dn’s,)  -  T,kX  tin*,] 

hJ  ’  b6  sina0()  -  VjkXsin*,) 

kj  "  b!  IT^kJ-kJdn1#,)  +  Tjk.k;**,! 

k4  -  bs  IV^kJ-kJsin1#,)  +  VjkX  sin«,l 

The  values  of  _An,  0n,  and  yn  will  now  be  written  in  terms  of  the  above  parameters: 


^1 

— 

hlhI 

s 

*>; 

ti 

S 

D’.- 

7L2 

= 

hlh; 

h 

s 

d1 

h 

m 

>;• 

^3 

* 

hlh; 

p* 

= 

d; 

h 

= 

o'; 

= 

h,h; 

Pa 

• 

D'l 

y< 

= 

o; 

- 

h2hj 

p5 

= 

d; 

y5 

s 

o'; 

^6 

= 

h2h*2 

Pa 

= 

d; 

y* 

« 

o'; 

-^7 

= 

h2h; 

Pi 

= 

d; 

?7 

s 

o'; 
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■4s 


4 

4 

.1 


Au 


*12  = 


*14  = 


*15  * 


The  basic  equation  for  <  AXA*  >  can  be  written  as 


<axa;>  = 


<1  +  Akltl  sin2«,) 


16  _ Pj_  (  '  -  g,P,  ♦  ' 

X,  (..♦  *„>(  (i+v„)  (l-v.)  T  1 


+  P„)e'up"*  (1+P„«t) 


(i-W 


I 
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Computer  Programs  for  Calculating  the  Backscattering  Coefficients 
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LIST  OF  SYMBOLS 


e(I) 

a{L) 


Position  vector 

Permittivity  of  the  random  medium 

Conductivity  of  the  random  medium 

Average  relative  dielectric  constant  of  the 
random  medium 

Average  conductivity  of  the  random  medium 


Standard  deviation  of  the  dielectric  fluc¬ 
tuations 


»»2 


Standard  deviation  of  the  conductivity  fluc¬ 
tuations 


L 

k3 

ei 

K 

X(i.i') 

G0U,r') 
J.  (-L) 

fi(X-l') 

I 


Mean  thickness  of  the  vegetation  layer 

Complex  propagation  constant  in  the  soil 

Angle  of  incidence 

Free  space  propagation  constant 

Infinite  space  dyadic  Green's  function 

Mean  wave  in  the  random  medium 

Infinite  space  scalar  Green's  function 

Scattered  electric  field  in  the  random 
medium 

Three-dimensional  Dirac  delta  function 
equal  to  6  (x  -  x')  6  (y  -  y')  5  (z  -  z’) 

Unit  dyadic 


Unit  vectors  in  x,  y,  and  z 
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Kronecker  Delta 


k  k  k 
ex’  «y’  •« 


Component  of  tie  effective  propagation 
constant 


2s(Jst>z) 


Value  of  k#t  for  horizontal  polarization 

Value  of  kei  for  vertical  polarization 

Two-dimensional  Fourier  transform  of  the 
scattered  electric  field  in  the  random  me¬ 
dium 


Ax(kx,ky),Ay(kx,ky),  Al(kx,  ky) 


Fourier  transform  of  the  amplitudes  of  the 
scattered  electric  field  in  air 


cot 
*x  »  *y  *  ** 


Illuminated  surface  area 

Correlation  distances  in  x,  y,  and  z,  re- 
pec  tively 


Backscatter  coefficient  for  volume  scattering 
and  for  the  case  of  horizontal  polarization 
transmit,  horizontal  polarization  receive 

Backscatter  coefficient  for  volume  scattering 
and  for  the  case  of  vertical  polarization 
transmit,  vertical  polarization  receive 

Backscatter  coefficient  for  a  randomly  rough 
surface  for  the  case  of  horizontal  polar¬ 
ization  transmit,  horizontal  polarization 
receive 

Final  backscatter  coefficient  result  that 
includes  both  volume  scattering  and  rough 
surface  scattering  ftr  the  case  of  horizontal 
polarization  transmit,  uorizontal  polariza¬ 
tion  receive 


o°  v  Final  backscatter  ooefTicient  result  that 

includes  both  volume  scattering  and  rough 
surface  scattering  for  the  case  of  vertical 
polarization  transmit,  vertical  polarization 
receive 

F  Fraction  of  water  by  weight  in  the  vege¬ 

tation 

Volume  of  vegetation  divided  by  the  total 
volume 

m  Standard  deviation  of  the  rough  surface  fluc¬ 

tuations  divided  by  the  correlation  distance 

f  Frequency 
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